Quantum Fractal Fluctuations 
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We numerically analyse quantum survival probability fluctuations in an open, classically chaotic 
system. In a quasi-classical regime, and in the presence of classical mixed phase space, such fluc- 
tuations are believed to exhibit a fractal pattern, on the grounds of semiclassical arguments. In 
contrast, we work in a classical regime of complete chaoticity, and in a deep quantum regime of 
strong localization. We provide evidence that fluctuations are still fractal, due to the slow, purely 
quantum algebraic decay in time produced by dynamical localization. Such findings considerably 
enlarge the scope of the existing theory. 
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The study of mesoscopic conductance fluctuations has 
gained important insights from the analysis of the un- 
derlying classical dynamics [Q. In particular, it has been 
found that some statistical properties of quantum fluctu- 
ations are semiclassically determined by the classical law 
which rules the decay in time of the survival probabil- 
ity inside the open chaotic system Q|. Building on this 
general relation, it was eventually predicted that, if the 
classical decay law is algebraic, due to the presence of 
residual stable islands in the classical phase space, then 
the quantum fluctuation pattern should display a frac- 
tal structure with a fractal dimension directly related to 
the classical decay exponent Such fractal structure 
of quantum fluctuations has been interpreted as a quan- 
tal manifestation of classical fractality associated with 
hierarchical structures of surviving stable islands. This 
theoretical prediction has received experimental support 
Q , and has been numerically confirmed on a dissipative 
model system, where coupling to continuum is simulated 
by complete absorption of that part of the wave packets 
which propagates outside the interaction region |^,^. 

In such a theoretical frame, quasi-classicality of the 
quantum dynamics appears essential in order that slow 
classical decay be quantally reproduced over sufficiently 
long time scales. In contrast, in a deep quantum regime, 
localization effects would become dominant, rapidly ef- 
facing any memory of the underlying classical dynamics. 

In this Letter we provide evidence that fractality of 
quantum fluctuations nevertheless survives even in such 
strongly non-classical situations. In the model systems 
presently to be described, the classical phase space has 
no significant stable islands, the classical motion is diffu- 
sive, and the classical decay is exponential; therefore, in 
the quantum weakly localized regime, one observes non- 
fractal, Lorentz-correlated conductance fluctuations [Q. 
In contrast, we work in a strongly localized regime, and 
observe fractal fluctuations, apparently unrelated to any 



classical phase-space structure. Decisive in this respect 
is the fact that the quantum decay, albeit highly non- 
classical, is algebraic due to purely quantum localization 
effects, as recently predicted |^,^. As a matter of fact, we 
find that the fractal dimension of quantum fluctuations 
is connected to the algebraic decay exponent in exactly 
the same way predicted by the semiclassical theory . 
It therefore appears that the quantum fluctuation pat- 
tern bears the same relation to the quantum probability 
decay, which was originally predicted by semiclassical ar- 
guments, even though such decay is not by any means 
quasi-classical. 

We consider the kicked rotator model with absorb- 
ing boundary conditions. Classical dynamics is described 
by the standard map: 



/ = / + A; sin 6*, 
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with absorption for / < and for / > A^. Classical dy- 
namics only depends on the scaling parameter K — kT. 
The corresponding quantum map is obtained by the sub- 
stitution I ^ fin = —ihd/dO. We denote -04 the state 
vector immediately before the t— th kick, T the kicking 
period, and we set fi = 1. Throughout this Letter, we 
will treat the kick counter t as a discrete time variable 
in units of T . The discrete-time evolution from time t to 
time t -I- 1 is described in the classical case by (|^) and in 
the quantum case by : 



(2) 



The state i/'t+i is obtained from i/'t by means of the uni- 
tary one-period propagator of the kicked rotor. The op- 
erator P then projects V't+i over states in a fixed interval 
of n values (0 < n < N), so it describes complete dele- 
tion (absorption) of the part of the wave packet which 
propagates outside the given interval. The parameter (p 
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can be interpreted as an Aharonov-Bohm flux through 
the ring parametrized by the coordinate 6. 

We consider the case K = kT = 7, k = 5, with absorp- 
tion for n < and n > N = 250. With such parameter 
values, the classical phase space has no significant island 
of stability We consider a statistical ensemble of 

orbits starting at t = with 7 = and randomly dis- 
tributed phases 6. The survival probability P{t) is the 
fraction of the ensemble which hasn't been absorbed at 
time t] equivalently, P(t) is the integrated distribution of 
exit times of orbits in the ensemble. 




FIG. 1. Classical (dashed line) and quantum (thick solid 
line) survival probability for the kicked rotator model with 
if = 7, fc = 5, = 0. The evolution starts at no = and 
the probability is absorbed for n < and n > 250. The thin 
solid line has slope one. Here and in the following figures the 
logarithms are decimal. 
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FIG. 2. Survival probability vs flux (f> at fixed time t — 10^ 
(above) and t — (below), with parameter values as in 
Fig.0. 

The classical decay of P{t) is well described by a dif- 
fusion process, practically insensitive to the chosen value 
of (f). As shown in Fig.^, after a diffusion time tjj w 
N^/D oc N'^/k'^ [D diffusion rate), the survival proba- 
bility exponentially decays with time: P{t) oc exp(— 7^), 



with 7 « l/^D ~ 2.2 X 10~^. In contrast, quantum dy- 
namics is strongly localized, as the localization length 
£ w fc^ is much less than the "sample size" N. Such 
strong localization causes a slow decay of the quantum 
survival probability P,f,{t) The latter is defined as the to- 
tal probability on states < 71 < 250 at time t, given 
that a,t t — the state was uq — 0. After the "Heisen- 
berg time" tn, it decays proportional to 1/t up to a time 
tmax ~ e-xp{2N/i) 1^, as shown in Figj^ (in this case, 

tmao; ~5 X 10«). 

Though the exponent 1 of algebraic quantum decay is 
independent of the value of (j), the quantum survival prob- 
ability P(ji{t) at a fixed time t is sensitively dependent 
on (p. The fluctuation pattern exhibited by the graph 
of PAt) versus (j) at fixed t (two examples are given in 
Fig.g) is precisely the object of our analysis, at values of 
t up to 10^, and for 10"' values of 4> in the interval [0, 0.1]. 
First we have computed autocorrelation functions: 



C((50) ^< P^{t)P^+s^{t) >^ . 



(3) 



Such correlations play a key role in the existing semi- 
classical theory of fractal fluctuations, which consists of 
three steps: 

(i) correlation functions are related via Fourier transform 
to the classical decay of the canonical variable conjugate 
to (j) H,^, notably, to the function Pci{Q) yielding the 
probability distribution of the angle Q accumulated by 
the rotor at the time of exit from the finite sample ( "accu- 
mulated" means multiples of 27r included). In particular, 
if Pc/(0) decays as at large 8, then 



C{S(j)) ~ C(0) - const|(5(/.|" 



(4) 



at small 5(j). 

(u) It is physically intuitive, and numerically confirmed 
0, that the decay exponent a coincides with the decay 
exponent of the classical survival probability, 
(iii) a signal-theoretic argument allows to conclude 
from (i) and (ii) that the fluctuation graph has a fractal 
dimension given by 



D = 2- 



(for a < 2), 



(5) 



where a is the time decay exponent of the classical sur- 
vival probability. Step (iii) actually invokes certain sta- 
tistical properties of fluctuations, which were assumed to 
be ensured by the chaotic nature of classical dynamics 

In our case, semiclassical arguments are invalid, be- 
cause of strong quantum localization. In addition, there 
is no classical decay exponent a, because classical de- 
cay is exponential. So the above theory doesn't apply; 
still, numerical results to be described below demonstrate 
that (i) , (ii) and (iii) are still legitimate in our highly non 
semiclassical case, too, provided one replaces the classical 
distribution functions P{t),Pci{Q) by the corresponding 
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quantum distributions P^{t),Pq{Q), and a by the expo- 
nent of the quantum algebraic decay. In order to compute 
Pq{Q) we perform a Ufting of our model (||), that is, we 
translate the model (||) of a kicked particle on a circle 
into the model of a kicked particle on the line, described 
by the coordinate 6. To simulate the corresponding dy- 
namics we restrict inside a large, yet finite box, with 
periodic boundary conditions. At every evolution step t 
we compute the function \ipt{Q) — ipt{'d)\'^, which yields 
the probability loss occurring at time t and position O 
due to absorption. Summing this quantity over t gives 
the total probability lost at position 6, that is, the prob- 
ability distribution of the angle accumulated at the exit 
time. The corresponding integrated distribution is shown 
in Fig.^ and displays a decay oc 8^^. The faster decay in 
the rightmost part of Figgis due to the finite total time 
and to the finite O— box used in the integration. Thus, at 
any given time t, the distribution of accumulated exit an- 
gles up to that time can be assumed to decay proportional 
to 6"^ in a range 6mm < Q < Om,ax,t, with &max,t in- 
creasing with t as long as t < tmax- One can estimate 
©mm as the angle accumulated until the time tu ~ 
This gives e„™ ~ /o"^ ^{D/T)TdT - i2/3)t%T. As- 
suming a maximal momentum ~ £ for the rotor started 
at riQ — 0, a, rough estimate for Qmax,t is provided by itT 
(©min ~ 2 X 10^ and Qmax.t ~ 2.1 X 10^ for the case of 
Fig.3). 




1 2 3 4 5 



log© 

FIG. 3. Quantum integrated distribution of the accumu- 
lated angle O at the exit time, computed over a total time 
t = 600, parameter values as in Fig.^. The straight line has 
slope one. 

Numerically computed correlations (j^) shown in the in- 
set of Figj^, at two different times, indicate that relation 
is still valid in our strongly localized case |jl^. We 
have computed the fractal dimensions of the correspond- 
ing graphs of P4,{t) versus (f) shown in Figj^ using the 
box-counting algorithm described in ^,|^ . The results are 
presented in Fig.^ and demonstrate that the graph has 
a fractal dimension 1.5 over a significant range of 6(j) 



scales, as predicted by equation (||). This range roughly 
lies in between scales 6(j)min,t, S4>max- The latter scale is 
inverse to ©mm. The former scale can be estimated by 
©maa; t (hcuce it dccrcascs as t increases), as long as it 
remains above a minimal i5(/)o ~ 27t/{N'^T) scale. This 
scale is imposed by the finiteness of the sample, and can 
be explained as follows. The evolution operator (^) has 
complex eigenvalues inside the unit circle. Varying cf) 
causes these eigenvalues to move, and their maximal shift 
under a change Scj) is estimated by TNScjj. Any time one 
moving eigenvalue comes close to the unit circle, a local 
maximum in the graph of P^{t) vs cj) appears. The num- 
ber of peaks produced by a single moving eigenvalue as </> 
changes from to 1 is at most ~ TN/{2tt). As there are 
N eigenvalues, the total number of peaks is ~ TN"^ /{2n), 
with an average (j>— spacing ^ ScjiQ = (27r)/(7V^r). As 
the fractal structure of the graph is due to the super- 
position of many tiny peaks originated in this way, no 
fractality can be expected below the S(po scale. 
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FIG. 4. Fractal analysis for the kicked rotator model, with 
parameter values as in Figjl|. The graph of P(f,{t) vs was 
covered by boxes of side Sep, and the largest excursion of 
P^{t) in each strip was recorded; summing over all strips 
and dividing the result by S(j>, we obtained N{S(j}). Data are 
shown for t — 10^ (circles) and t — lO** (diamonds). The 
straight lines correspond to the fractal dimension D = 1.5 
expected for a = 1. The inset shows the correlation function 
C{5<j)) at the same times. The straight lines correspond to 
1 - C(5(^)/C(0) oc |(5(/>|", with Q = 1. 

Further confirmation for the validity of the above illus- 
trated scenario was obtained by analyzing the following 
variant of the basic model (H): 

which differs from (||) because the kicking potential has 
now a discontinuous derivative. Localization is in this 
case algebraic, as eigenfunctions decay like l/n?. The ar- 
gument used to predict the decay of the survival prob- 
ability in the case of exponential localization [B yields an 
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asymptotic decay oc t~^/^. The numerical analysis illus- 
trated above for the case of exponential localization was 
replicated for this model, too, yielding the results shown 
in Fig.||. Fractal analysis of the fluctuation graph yields 
good agreement over a broad range with the theoretical 
D value of 13/8 predicted using equation (||) and the 
decay exponent a = 3/4. 
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FIG. 5. Same as in Fig.|l| (upper inset) and Fig.^ but for 
the discontinuous map (^. The straight lines correspond to 
a power law decay of the quantum survival probability with 
a — 3/4 (upper inset), a correlation function C{5<j>) oc |(5(^|" 
(lower inset) and a fractal dimension D — 2 — a/2 — 13/8. 
The a;-axis in the lower inset is the same as in the main figure. 

In conclusion, we have provided numerical evidence for 
fractal fluctuations of the quantum survival probability 
in a classically chaotic system in the regime of strong lo- 
calization, and in the absence of significant classical crit- 
ical structures. Whereas the original theory |^ predicted 
fractal conductance fluctuations, our present results are 
about fluctuations of a different quantity, namely sur- 
vival probability. As both quantities ultimately reflect 
the fluctuations of the scattering matrix, their fluctua- 
tions should be similar in nature. As a matter of fact, 
the very same semiclassical arguments used to predict 
conductance fluctuations could be adapted to the sur- 
vival probability, too The main difference is that a 
scan of (t) at finite however large t cannot achieve the 
0— resolution exhibited by conductance at fixed (quasi- 
)energy. On account of the effect of increasing t reported 
earlier in this Letter, such coarsening is unlikely to "frac- 
talize" otherwise non- fractal fiuctuation patterns. In any 
case, survival probability is a meaningful quantity in a 
broader class of problems than electronic transport in 
semiconductor structures. The scenario we have ana- 
lyzed in this Letter can also be exported to other realistic 
problems, where 'conductance' is instead a problematic 
concept. As an example we quote microwave ionization 
of Rydberg atoms |^ . This may open a way for experi- 
mental observation of localization-induced fractal fluctu- 



ations in atomic physics. 

Our results show that quantum fluctuations may be 
fractal even in situations where the existing semiclassical 
theory does not apply. This does not command a reinter- 
pretation of existing experimental data, which were ob- 
tained in situations far from the strongly localized regime 
considered in this Letter. Nevertheless it signals that the 
current understanding of the deep quantum mechanisms 
responsible for fractal fluctuations is still far from com- 
plete. 
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